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Abstract 

We extend the calculation by Turfus and Schubert (2016) of analytic pricing of CoCo bonds under a 
Hull- White short rate model for equity conversion intensity to a beta blend short rate model which 
encompasses both Hull-White (normal) and Black-Karasinski (lognormal) models. The asymptotic 
expansion assumes that the volatility of the conversion intensity is small. Results are calculated up 
to second order accuracy in the volatility. This is considered adequate for most practical purposes. 


1 Introduction 

A hybrid credit-equity model was proposed by Turfus and Schubert (2016) for the pricing of contingent 
convertible bonds (CoCos) with reduced form modelling of the occurrence of conversion events. The reader 
is referred to that paper and references therein for details about the defining characteristics of CoCo bonds 
and the various pricing approaches which have been proposed in the literature. The evolution of the equity 
conversion intensity was assumed to be governed by a diffusive process, and the price of the equity to which 
the bond converts by a jump-diffusion process, with a downward jump of a fixed relative amount k (< 0) 
occurring at the time r of the conversion event. Thus the equity payoff was taken as M( 1 + k)S~, with M 
the contractual conversion rate. The diffusive processes are further assumed to be (negatively) correlated. 

The main achievement of the paper was to derive, by means of a perturbation expansion approach, 
approximate analytic solutions for the hybrid model under the assumption that the credit /conversion process 
was specifid by the (normal volatility) short rate model proposed by Hull and White (1990). Here we make 
two extensions to this work. Firstly, we embed the Hull- White model within a more general framework that 
includes in particular the (lognormal volatility) model of Black and Karasinski (1991). The latter has the 
advantage that it ensures intensities remain positive, which the Hull- White model does not. Secondly, we 
extend the first order approximation of Turfus and Schubert (2016) to include second order terms. These 
turn out, as we shall see, to be more important in the Black-Karasinski (lognormal volatility) case than they 
were for the Hull- White case. 

We also note here that the model we propose here has many similarities than one of two models proposed 
by Cheridito and Xu (2015). They consider a structural model which seeks to represent the capital ratio 
of the CoCo-issuing bank directly, with conversion occurring when it drops below a trigger level. They also 
propose a reduced form model with many similarities to the one we propose here. Their models are solved by 
application of finite differences to the governing PDEs. Three significant differences between their reduced 
form model and the one we set out below are as follows: 

• While we ensure positive intensities by using a Black-Karasinski model, they do so by using instead 
the square root model of Cox, Ingersoll and Ross (1985). 

• While they take the payoff in the event of a conversion event at time r to be M(l+k)S~ unconditionally 
(in our notation), we suppose this amount to be capped, at least by the bond notional, on the basis 
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that it is highly improbably the capital ratios of a bank would trigger a conversion event when the 
equity conversion price exceeds the value of the convertible bond. 

• Finally, Cheridito and Xu (2015) take default events explicitly into account insofar as they consider 
there to be a small probability that, when conversion is triggered, this does not save the bank from 
default, whence, rather than the stock falling to (1 + k)S~ , ft falls to zero. 

The first of these is a residual difference which probably ought not to have too much impact on the conclusions 
reached since, as we shall see, the conversion intensity model only impacts through the integrated term 
variance and covariance of the intensity, which properties are likely to be similar under parallel calibrations 
of different models to the same market data. Our use of a possible cap makes our model more general than 
that of Cheridito and Xu (2015). We shall return to the final difference below but, suffice it to say, the 
incorporation of the impact of default events turns out to require only a trivial modification to our model. 
On this basis we suggest that all the most important features of the model of Cheridito and Xu (2015) are 
included in our model, with the advantage that analytic formulae are available to compute prices and risk 
sensitivities to a good degree of accuracy. 

The structure of the paper is as follows. The founding assumptions of the mathematical model are set out 
in §2 below. The perturbation analysis leading to the calculation of the equity recovery value of the CoCo 
bond is set out in §3, the main results of which are Theorem 3.1 and Proposition 3.1. The accuracy of the 
asymptotic results and the importance of second order terms are considered briefly in §4. Consideration is 
given in §5 to how the impact of explicit default events can be taken into account in the manner of Cheridito 
and Xu (2015). 

2 Mathematical Model 

We consider an equity conversion process driven by a (stochastic) conversion intensity A t , the dynamics of 
which are governed by a beta blend short rate model. The method follows closely the similar calculation 
performed by Turfus and Schubert (2016) using a Hull- White short rate model. The equity price is taken to 
be given by a jump-diffusion process, with a downward jump of a fixed relative amount occurring at the time 
of a conversion event. The diffusive processes are potentially correlated. It having been observed by Turfus 
and Schubert (2016) that stochastic interest rates have little impact on pricing (at least in the non callable 
case), we assume from the outset that rates are deterministic. 

We shall find it convenient to work not with the conversion intensity X t directly, but with a reduced 
variable y t satisfying the following canonical Ornstein-Uhlenbeck process: 

dy t = -ay t dt + <J\{t) dW (1) 

where a is a positive constant specifying the mean reversion rate, 1 cr\(t) is a bounded positive L 2 function 
and dWf is a Brownian motion for t £ D m := [to, min{r, T m }} where r is the time of the first relevant CoCo 
bond conversion event and T m is the longest maturity date for which we wish our model to be calibrated. 
This equation is well known to have a unique strong solution subject to these assumptions. The variable y t 
is taken to be related to the conversion intensity X t by 

(1 -0)\ t + /3\(t) - A* (t) exp ( ( ) (2) 

with /3 £ [0, 1) assumed and A (t) > 0 the instantaneous forward rate determined at the initial time to with 
reference to market data and a suitable interpolation methodology. Here A*(t) is chosen so as to satisfy the 
no-arbitrage condition that 

7-t T — If A s dsl — A (s)ds /0 , 

E e o = e Jt o w (3) 

1 It is a straightforward matter to loosen this assumption and specify instead that a(t ) be a bounded positive L 1 function, 

whereupon the analysis below goes through effectively replacing a(v — u) throughout by f” a(s)ds. 
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under the martingale (money market) measure for to < t < T m . It is not difficult to show that Eq. (2) is 
consistent with Eq. (1) and furthermore that the local volatility of X t is 

(l-f))A ,+0X(t) 

— w — n(t) - 

This can be thought of as a displaced diffusion-type interpolation between the Hull- White (normal) and 
Black-Karasinski (lognormal) models. For the equity process we propose: 

7 Q 

— ^ = (r(f) - q(t.) - k\ t ) dt + a s {t) dW f + kdn t , (4) 

where S t is the equity price, q(t) its expected (continuous) dividend rate, f{t ) the instantaneous forward rate 
of interest, dW% is a Brownian motion satisfying 


corr (W t \W?)=p xs 


(5) 


and nt a Cox process with intensity A t giving rise to an equity price jump of size k with — 1 < k < 0, 
contingent on an equity conversion event. 

To proceed, we follow Turfus and Schubert (2016) and seek to price the survival-contingent flows 
separately from the conversion-contingent flows, expressing the PV of a CoCo bond with notional N as 


V(t) = Nf SUIV (t) + Nf conv (t). 


( 6 ) 


where / Sllrv (t) expresses the part of the PV generated by survival-contingent cash flows, i.e. the coupon 
payments and notional repayment, and f CO nv(t) expresses the contribution to the PV resulting from cash 
flows contingent on a conversion event at time r < T, the bond maturity. We focus here on the calculation 
of the latter term and refer the reader to Turfus and Schubert (2016) for that of the former. 

We propose as our model for the payoff of the CoCo bond upon conversion at t = t that this is given 
by min{A/S'“(l + fc), A'}, with M the number of shares issued per unit notional on conversion, and K < 1 
specifying a cap on the value of the equity payoff. 

Following Turfus (2016), we seek an approximate solution f CO nv(t) under a “weak volatility” assumption. 
To this end we rescale both y t and <7\{t) by an asymptotic parameter e defined by 


e 2 := 


^(*) 


dt 


a(T m ~t 0 )J tn \X(t)\W 

which we take to be small. Thus we define new scaled variables y t and <J\{t) by 


Vt ■= e 1 yt, 
a y (t) := e~ 1 a\(t), 


both taken to be 0(1) as e — > 0. 

We further define an equity term variance 


Is(t. i,t 2 ) := / cr| (u)du, 


a new characteristic stochastic equity price coordinate Xt such that 

S t = F(t)e Xt ~^ Is(t 0’*) 


(7) 

( 8 ) 


(9) 


( 10 ) 


where _ 

F{t) ■= S t 


(ID 
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a payoff function 

M 0 (x,t) := M{l + k)F(t)e x ~^ Is{to ’ 1) 

and a risky discount function 

B(t u t 2 ) := e~ Sl?<s(e)+\{B ) )ds . (12) 

We further suppose that the process f conv (t) can be expressed in terms of the new (stochastic) variables a : t 
and y t as h(x t ,y t ,t). We are interested in calculating f conv (t o) = h(0,0,t o ). 

Applying the well-known Feynman-Kac method to Eqs. (1), (4) and (5), we infer that the function 
h(x, y , t) is governed by the following backward diffusion equation: 


C[h(x,y,t)\ 


(A (t) + AA (y, t)) min{M 0 (a;, t), K} + AA (y, t) 



(13) 


for t £ D m with final condition that h(x,y,T) = 0, where £[.] is a standard forced diffusion operator given 

by 2 2 

41 [ =Ft + 1 2 + ^ (t) |?) “ (f(t) + J(t)) (141 

and 

AA to a):= r ^(v W e x p(Agk)_X ( t) ), (15) 

so that A t = A (t) + AA (y t ,t) 


3 Perturbation Analysis 

In the absence of exact closed form solutions to Eq. (13), we seek a solution as a perturbation expansion in 
e. To that end we write _ 

X * { llp {t) = eA * W + ^2 W + + 0(e 4 ) 

where we know from Turfus (2016) that the no-arbitrage condition Eq. (3) requires us to choose A *(t) = 
Ag (t) =0 and 

A = f Jt q,u) du- - ® J y (to,t). (16) 

| A(t ) Jto |A(«)|^ 2 7 |A(<)I 2/3 

with t 

Iy(ti,t 2 ) := f 2 e- 2a ^- u) a 2 y (u)du. (17) 


By expanding the exponential in Eq. (15) as a power series about e = 0 and gathering together like powers, 
we can also write 


AA (y, t) = eAi (y, t) + e 2 X 2 (y, t) + 0(e 3 ) 

(18) 

where 


V 

Ti- 

ll 

(19) 

= ^(t)+ 4 (l 0) y 2 . 

(20) 

Note that by construction we have A* ( t ) = Aj(0, t), i = 1, 2, . . . 
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To solve Eq. (13) in it asymptotic representation we pose a perturbation expansion 


h(x,y,t) = h 0 (x,t) + ehi(x,y,t) + e 2 h 2 (x,y,t ) + 0(e 3 ). 


( 21 ) 


Substituting the above expansions into Eq. (13), solving for successive power of e, setting t = to and reverting 
to unsealed notation, we obtain the following first order accurate result: 

Theorem 3.1 The value of the equity recovery payoff on a CoCo bond can be estimated under our modelling 
assumptions as follows: 


fconv(t 0 ) = J B(t 0 ,v)X(v)(M(l + &).F(t;).M[lV'(-di(a:,to,v))] ^ + KN(d 2 (0, t 0 , v))) dv + 0(e 2 ) (22) 


where the operator M is defined by 


A4[f(x)] \=(\ + _ f j ^ ^ ^( u ) , A _[_ fc _[_ /(%) 

V |A(n)|/5 J t0 Ry °’ |A(^)|/3 V dxJJ ’ 

and 

Iit(ti,t 2 ) := pxs f e~ a{t2 ~ u) cr\(u)<Js(ii)du, 

Jt i 

In M 0 (a:, u) — In if 

d 2 {x,t,v) := , 

V I s(t,v) 

di(x,t,v) := d 2 (x,t,v) + y/l s (t,v), 

1 r x 

N(x) := / exp (— Am 2 ) du. 

v 27T ,/ — oo 

Proof. For the proof of Theorem 3.1, see Appendix A. 


(23) 

(24) 

(25) 

(26) 
(27) 

□ 


We make the observation that the form of Eq. (22) is akin to a CVA calculation, giving the PV of a 
payment at default of an equity underlying capped at a collateral level K. Ignoring the second two terms 
on the r.h.s. of Eq. (23) and setting M = 1 gives the value ignoring wrong-way risk and jump risk. The 
neglected terms correct for the latter effects. It is worth observing that, although the wrong-way risk is 
initially (for short times to maturity) driven primarily from the second term in direct proportion to the 
covariance Ir(-), as time to maturity is increased and/or the default /conversion intensity increases, the 
impact of the covariance is diminished by the impact of the remaining third term. 

Although the result of Theorem 3.1 gives the equity recovery value to first order in e, it is at times 
necessary to use a more accurate approximation, particularly for the Black-Karasinski case where <J\{t) will 
not in general be particularly small. This can become an issue when the time to maturity of the bond exceeds 
about 10 years. Errors can exceed about 50 bp of notional rendering the first order approximation less than 
adequate (see §4 below). To that end, the process described in Appendix A can be extended to allow second 
order terms to be calculated. The following result is obtained. 


Proposition 3.1 A second order accurate estimate for f CO nv(t o) is obtained by adding to the contribution 
of Eq. (22) the following adjustment: 

Afconv(to) '■= f B(to,v)X(v)M(l + k)F(v)ip(0,to,v)dv (28) 

Jto 

Details of the derivation of this expression and definition of ip(.) are provided in Appendix B. As can be 
seen from the form of a quadratic dependence on the correlation level p\g is introduced at this order of 
approximation through Ir(-), and an explicit (quadratic) dependence on the conversion intensity volatility 
level, through Note that the contribution from the latter vanishes when the jump parameter k — > 0. 
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We observe that the Hull- White case (/3 = 1) was excluded under our initial assumptions as a singular 
limit of Eq. (2). However it will be observed that /3 can conveniently be set to unity in the above expressions 
f° r f com (to) an d A/ C0nv (fo)- It can be confirmed that the expressions deduced are consistent with those 
obtained by Turfus and Schubert (2016), leading us to conclude that the results presented herein are indeed 
valid for all /3 £ [0,1]. 

4 Results and comparative analysis 

Computations were made for the important Black-Karasinski case of the equity recovery value of a range of 
CoCo bonds under various market assumptions based on the zeroth, first and second order expansions. We 
chose values of A = 6%, as = 30%, a\ = 50%, a = 35%, p\s = —0.6 and k = —0.8. The notional N was 
taken to be 100; the initial share price S to , the share issue rate M and the cap I\ were all taken to be 1. 
The results are compared in Figs. 1 and 2, with the time to maturity in Fig. 2 taken to be 15 years. 



-0th order 
-1st order 
-2nd order 


Figure 1: 


Asymptotic estimates of equity recovery value for various CoCo bond maturity dates. 
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1st order 
2nd order 


Figure 2: Asymptotic estimates of equity recovery value with various assumed conversion intensity- 
equity correlation levels. 

It is observed that, while the second order corrections for maturities of 10 years or below were less 
than 20 bp of notional, this quickly increased to 50-60 bp for maturities of 15y or more for typical market 
circumstances. This is larger by a factor of about two than the size of the errors observed by Turfus and 
Schubert (2016) for an assumed normal conversion intensity volatility. It is therefore suggested that, for 
longer maturity bonds, the second order expansion should be used in preference to the first order. However, 
given the large decrease in size of the second order correction compared to the first order, there is good 
reason to believe it justifiable to ignore the contribution due to third and higher order terms and to suppose 
our second order expansion to be adequate for most practical purposes. 

5 Extensions of the Main Result 

Cheridito and Xu (2015) in their reduced form model propose an additional feature not considered in 
the above model, namely the possibility of a default being triggered by the conversion event. They take 
the probability of this happening (with zero resultant conversion payment) to be a (<C 1), whence the 
probability that a payoff of M (1 + k)S~ is received is 1 — a. But this feature is easily incorporated into our 
framework: instead of a deterministic jump of fc, we face a binary choice between a jump of k and a jump of 
— 1. Consequently, rather then the (uncapped) payoff being M(l + k)S~ , we must consider it to be effectively 
given by (1 — a)M(l + k)S~. Having made this change, to ensure that the cap is applied correctly to the 
effective payoff, we should set it to (1 — a)K rather than K. Finally, rather than the jump-compensating 
drift of — k\t in Eq. 4, we need to use a drift of — ((1 — a) x k + a x — l)At = — (k — a(l + k))Xt- In conclusion 
our model only needs to be modified by changing 

k —> k — a( 1 + k) 

K (1 - a)K (29) 

to incorporate the impact of default at conversion with probability a. With this change, all our main results 
can be used as presented. 

We add in passing that it was observed by Turfus and Schubert (2016) that it was necessary to choose 
rather higher values of k than might have been expected a priori to achieve a consistent calibration of the 
model to market prices of CoCo bonds. This was interpreted by suggesting that the large jump apparently 
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being priced in was not just from the impact of the conversion event itself but from other factors impacting 
the stock price at around the same time as the conversion event. The occurrence of a default immediately in 
the wake of a conversion event is a good example of such. On this basis we propose that factoring in a small 
probability of “default at conversion” is probably an appropriate measure to take in looking to calibrate the 
above model. 

We make the further observation that Cheridito and Xu (2015) propose that, if there is no immediate 
default upon conversion, there remains a risk of a subsequent default, with the relevant intensity given by 
fit = (3\t for some constant 0 < /3 < 1. The value of this constant does not impact on the pricing but is 
used by the authors in their calibration process, insofar as they look to have the model price CDS rates in 
a manner consistent with this double-event model. By assuming a and /3 to be known a priori, they are 
able to use the CDS calibration process to infer the (otherwise) unknown conversion intensity A*. While 
this is superficially attractive, it is arguably only exchanging one intractable problem for another. For, if we 
suppose the conversion intensity to be intrinsically unknown, we are not in a position to infer its relationship 
to the default intensity, other than that it will be greater. The approach of Cheridito and Xu (2015) is 
essentially to suppose the following: 


a\ t , t <t, 
PX t , t > T , 


in other words a piecewise linear relationship. While this ansatz would appear to be capable of yielding 
information about A t from the better understood y t , confidence in any conclusions drawn is only merited 
to the degree that a and p are known. Of course, under the approach of Cheridito and Xu (2015), the 
assignment of both values is necessarily an ad hoc process. In the end, many more equally plausible calibration 
approaches could be used, such as inferring directly from CDS and postulating that A t = Ht + s(t), for 
some deterministic function s(t), chosen to attain consistency with CoCo market prices. 


A Proof of Theorem 3.1 


We derive Eq. (22), solving Eq. (13) in the standard manner by successive levels of approximation. At zeroth 
order we must solve 

JZ[ho(x,t)] = — X(t) min{Mo(a:, t), K} 

for t. £ D m with final condition that h${x,T) = 0. This can be achieved by means of the following readily 
obtainable Green’s function for the canonical diffusion operator £[.]: 

G(x, y, t; £, r/, v) = B(t, v)H(v - t)cj)(x - £, ye~ a(v ~ t] -r];R(t,v)) (30) 


where B(t,v) is the risky discount factor defined in Eq. (12) above, H(.) is the Heaviside step function and 
(j)(xi, X 2 ] R{t, v)) is a bivariate gaussian probability distribution function with mean 0 and covariance matrix 


R(t, v) 


Is(t,v) I p (t,v) \ 
I p (t,v) Iy(t,v) ) 


(31) 


where 


= Pxs [ e o(t2 u) <j y (u)a s (u)du, 

Jt i 


and Is (t\ , £ 2 ) and I y (t\,t2) are defined in Eqs. (9) and (17), respectively. Making use of Eq. (30) we obtain 


/»oo />oo 


h 0 (x,t) = 


' t J — OO J — OO 


G(x, y, t ; £, 77 , v)X(v) min{Mo(£, v), K} c?£ dy dv 
X(v)B(t,v)(e*-i Is ( to ’ t) M(l + k)F(v)N(-d 1 (x,t,v)) + KN(d 2 {x,t,v)j^dv (32) 
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where di(.), d 2 (-) and N(.) are defined in Eqs. (25 )- (27) above. Continuing to first order we find 


C[hi{x,y,t )] = Ai (y,t) min{M 0 (x, t),K} + h 0 + 

for t £ D m with final condition that hi(x,y,T) = 0. The solution can be written 

hi(x,y,t) = h lfi (x,t) + yhi,i(x,t), 

where, applying our Green’s function as previously and integrating over 77 , we find 


v,N '(zmui) Ie(t ' v) iX 


\Hv)\p 


^min{M 0 (£, v), K} - ^1 + k — j h 0 (£, v)j d£ dv 

= e x ~i Is ( t °’ t) M(l + fc) ^ F(v)B(t,v)I p (t,v)j^^N(-d 1 (x,t,v))dv 

- e ^-|^(to,t) M ( 1 + J F(v)B(t, v)X{v) 

f I p (t, u ) ( 1 + k + k-^- ] N(-di(x + I s (t, u),t, v)) du dv. (33) 

Jt |A(u)|p V ox ) 


Likewise we find 


hi t i(x,t) = J B(t,v)N'^~ 


,( ^~ x \ -afv-t) Mv) 


/ Is{t,v) 


|A(v)|^ 


^min{M 0 (£, v),K} - ^1 + fc — J /i 0 (£, v)j d£ dv 
= [ B(t,v)e- aiv ~ t) 3 V \ (e x -^ Is{to ^M{l + ^F^Ni-d^x^^)) + KN{d 2 {x,t,v))) dv 

Jt | A (^)!^ 3 V > 

- f T Bit , v)X(v) [ V -^ 7 - + fc) 2 F( V )lV(-d 1 (x + I s (t, u), t, v)) 

Jt Jt |A(w)r v 

+ KN(d 2 (x, t, v))'jdii dv. (34) 

Substituting the above expressions into Eq. (21), setting t = to and reverting to unsealed notation, we obtain 
Eq. (22). 

This concludes the proof of the theorem. □ 


B Derivation of Second Order Terms 

For notational convenience we define at this stage the following differential operators: 

A, 


D\ := 1 + J-, 

OX 

d 

D 2 •— 1 + , 

OX 

D 3 •— 1 k k- 
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and an unsealed equity conversion variance 

■= e 2 I y {ti,t 2 ). 

Continuing in the same vein as Appendix A we find that the second order contribution h 2 (x,y,t) satisfies 
the following: 

C[h 2 (x,y,t)\ = \ 2 (y, t) (— min{M 0 (x,t), K} + D 2 h 0 (x,t)) + Ai(y, t) (p 2 h 3fi (x,t) + yD 2 h 1A (x,t)j , (35) 

with final condition h 2 (x, y, T) = 0. From the form of Eq. (35), we are led to suppose a solution of the form 

h 2 {x , y , t) = /i 2 , 0 (x, t) + yh 2fl (x , t) + y 2 h 2 , 2 (x, t), (36) 

Applying our Green’s function Eq. 30 in the standard way, carrying out the required integrations and 
simplifying, we find in particular that 

rT 


2 h 2t o(x,t) = e x 2 J s( t o.*)jvf(i —|— fc) J F(v)B(t,v)\(v)ip(x,t,v)dv 


(37) 


where, reverting for convenience to unsealed notation, we have 

t2u^ t.u rv J(u 


ip(x,t,v) := N{-di{x,t,v)) - 


- / e 


-±(l-P)J t 

r Xu 


|A(«)|^ 


\Xu )\ 20 

Xu) 

\\( u )\ 2 0 
lR(t,u) I 


I R (t,u) du D 3 N(-di(x,t,v)) 
(kl\(t,u) + I%(t,u) ((1 + k)D 2 - D%) ^N(-di(x,t,v)), 

(1 + k)Ift(t,u) du (. D 2 - Dl)N(-di(x,t,v)) 


\\(v)\ 2 « |A(u)|/3 J t |A(u)|* 

— a(v—u ) 


A (w) 


+(1+k) rs 

Jt |A(u)|£ 


|A(w)|^ 

p —a(u—w ) X(lV 


I R (t , w) dw du D 2 DiN(—d\(x, t, v)) 


|AHI* 


( 'kI\(t,w)D 3 + Ix(t,w) ((1 + k)D 3 Dj - D 2 Dl ) ^N(—di(x,t,v)) dwdu. 


(38) 


The final 0(e 2 ) contribution to / C0nv (io) is then given by e 2 /i 2 (0, 0, to)- Clearly the y-dependent terms in 
Eq. (36) do not contribute (whence we do not produce them here). Substituting in from Eqs. (37)-(38) we 
obtain the result in Eq. (28). 
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